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Abstract. We obtain the asymptotic distribution of eigenvalues of real 
symmetric tridiagonal matrices as their dimension increases to infinity and 
whose diagonal and off-diagonal elements asymptotically change with the 
index n as Jnt+i nt+i ~ aiV?(n), Jnt+i nt+i+i ~ h^{n), i = 0, 1, . . . , t - 1, where 
and hi are finite, and belongs to a certain class of nondecreasing 

functions. 
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A very versatile method for calculation of physical properties of quantum 
systems and references in 0-0) uses the Lanczos procedure to reduce 

the matrix of the Hamiltonian or some other relevant operator to the tridi- 
agonal form. After this reduction, various properties of tridiagonal (Jacobi) 
matrices provide much help in solving the problem at hand. An important 
problem arising in connection with this method is to find the asymptotic dis- 
tribution of eigenvalues of Jacobi matrices (when the dimension of the matrix 
increases to infinity) if the asymptotics of the matrix elements is known. The 
purpose of this note is to present a solution to this problem in a particular 
case. 

Let J be a symmetric Jacobi matrix with real matrix elements (that is 
Jik = Jki, i,k = 0,1,2,..., Jik = if |2 - A;| > 1, Ja+i ^ ). Suppose 
that limn-^^Jnt+int+i/f{n) = ai and limn^oc Jnt+int+i+i/ f{n) = k, i = 
0, 1, . . . , t — 1, where Oj, 6j < oo are not all simultaneously zero; t is a positive 
integer; the function (f{x) : R"'' is nondecreasing and satisfies the 

condition lim^^oo V'(^ + x)/ip{n) = 1 for any real x. Note that for (p{n) = 1 
the matrix J is asymptotically periodic. 

Furthermore, let J{n) = \ \Jik\\^^k^o denote the nt x nt truncated J. Our 
aim is to calculate the eigenvalue density p{z) of J{n)/ip{n) as n — > oo. It is 
defined by the requirement that p{z)dz be in the limit n ^ oo the relative 
number of eigenvalues of J{n)/(p{n) in the interval dz. 

We now impose a further limitation on the form of ^p{n): suppose that 
there exists a continuous function defined by g{uj) = du/du, where u^u) = 
limn^oo fi'nu) / (f{n) for u G [0, 1]. Note that for f{n) = rC , 7 > 0, we have 
g{u) = c.-i+V7/^. 

The ideas of these definitions come from the works P,^, where the case 
t = 1 was considered, and the moments of the eigenvalue density p{z) were 
calculated. The explicit expression for p{z) was given in for t = 1, a = 0. It 
is called Nevai-Ullman distribution. Using a different approach. Van Assche 
has obtained for ip{n) = n'^, and arbitrary a, b the density p{z) for t = 1 []TU| 



and (implicitly) for t = 2 [|TT[]. We will now obtain p(z) for any t and a, b by 
generalizing the method of |^,|^. 

Before that, let us recall some facts about periodic symmetric Jacobi 
matrices. Let L be such a matrix defined by the formulas Lnt+i nt+i = ca and 
Lnt+i nt+i+1 = bi, i = 0,1, ... ,t - 1, n = 0,1, .. .. The same a^, k appear in 
our definition of J. Define two systems of polynomials Pkix) and qk{x) by 
the initial conditions po = 1, pi = {x — ao)/6o; Qo = 0, gi = —bt-i/bo and the 
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same (for yk = Pk or yk = qk) recurrence relation: 

yk+i = {{x - ak)yk - bk-iyk-i) /h, k 



1,2, 



Let S{x) = pt{x)+qt-i{x). The spectrum of L consists of exactly t bands: the 
image of [—2, 2] under the inverse of the transform S{x) = A. The boundaries 



of the bands (denote them /ij 



t, 1^0 



-OO < /Ui < Z/i < /i2 < 



1/2 < ■ ■ ■ < Ut < IJ't+i = oo) are solutions to the equations S{x) = ±2. Thus, 
the spectrum of L is spec(L) = U-^]^[/ij, z/j]. The eigenvalue density of L is 
given by the formula po{x) 
and po{x) = otherwise. 

' -1 



\dS{x)/dx\/{t7i^4- S{xy), for x e spec(L), 



Note that for t = 1 we have po{x) 



[X 



ao) 



Theorem Let f{z,uj) = g{uj)po{z/uj)/uj. Then, under the conditions 
specified above, the eigenvalue density of J{n)/ip{n) as n ^ oo is the follow- 
ing: 

1) E [pk, t^k] for some k (1 < k < t). 



Iz/ 



f{z,u;)duj, 



_|_ y^t r^/Pj 
z/ui ' /-^j=i+l Jz/uj 

the preceding expression with i 

T.u~^j:',:;]f{z,u)du, 

f{z,uj)du, 



rl _|_ v^i— 1 r^/'^j 



the preceding expression with i 
^^'"'^ f{z,uj)du, 



z/nj 



z G [pi.Pi], Pi>0 

k, z e [0, z/fc] 

z G {ui,pi+i), Ui>0 

z e [pi, Vi], z/j < 

k, z e [pk, 0] 

z G {vi,pi+i), pi+i < 



2) E {vk, Pk+i) for some k (0 < k < t). The expression for p{z) is the 
same as above (two formulas for z G [0, z/^] and z G [pk, 0] should be dropped) 
except in the interval {vk^ Pk+i) where 

^■z/l^Lj' 



/-ij=k+l Jz/uj 



Z/Vj 

zliij 



f{z,uj)duj, 
f{z,u)duj, 



z G (0,/ifc+i) 
z G (z/fc,0). 



Sketch of the proof. Some basic ideas of the proof we borrow from 
where the case t = 1 was considered. Let us denote Km = I^oo^'^ Poi^)dx, 
where M is a nonnegative integer, and verify first the following identity: 

M \ 

Cfct+i, Gkt+i = Km J (1) 



1 

lim - V" 

k^oo t -0 



J 
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where {ej}°2,o is an orthonormal basis in which the above defined matrices 
of operators J and L are written. Indeed, it is obvious that the Lh.s. of this 
identity is equal to 

1 

hm -^(L^'^efct+ijefct+i), 



fc^OO t ■ n 



while 



t— 1 nt—1 

lim - Y^{L'^'' Ckt+u ekt+i) = 1™ ^ H (^^e^, ej) = Km- (2) 
From (|l]) it follows that 

1 n— 1 t— 1 n— 1 

- Y: E('^*'eH+., eH+.) = E ^(^)''(^M + e,), (3) 

k=Oi=0 k=0 

where — > as A; — > oo. We have: 



?i— >oo 



U,„ = 1, Un, = 0. (4) 



Hence, can be rewritten as follows: 

1 nt—1 „i 

- E (J^'ek, ek) = {Km + 5„)n / v{xnf'dx, (5) 
^ fc=o 

where 5„ ^ as n oo. 

To proceed further, we will need to estimate the sum 

nt—1 nt 

Zn{M)=Y.{J^'eu.eu)-Y.<n. (6) 

k=a k=l 

where Xk,^ k = l,...,nt are the eigenvalues of J(n). Note that Z„(M) 
is the difference between traces of two matrices, one of which is obtained 
by first taking the power J*^ and then truncating the result; and the other 
one, by first truncating J and then taking the power J{n)^ . It is easy 
to calculate that for n large enough Z„(0) = Z„(l) = 0, Z„(2) = Jntnt-ii 
^n(3) = Jit nt-ii.Jnt nt + Jnt-1 ni-i), etc. In general, we have for all sufficiently 
large n: 

|Z„(M)| </(M) max |J,,r, (7) 

\)<i,k<n+M 
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where f{M) does not depend on n. In view of this estimate, (^) leads to the 
identity: 

.n't Tm, = (8) 
nt Jq (^{xnY^dx 

Changing the variables Xk,n = ^in)zk,n and taking the limit, we have: 

^ rl pi poo 

rriM = lim — ^ = / uj^\x)dxK m = / uj'^^ g{Lu)duj / x^^ po{x)dx, 
M = 0,1,..., 

(9) 

where u{x) and g{x) are defined as above. 

Changing the variables z = ujx, ip = u in the double integral in (|^), we 

get: 

/oo 
z^^p{z)dz, (10) 
-oo 

where is defined in cases 1 and 2 as in Theorem. Since the actual support 
of integration in (|TU|) is finite, the function p{z) is defined uniquely by rriM, 
M = 0, 1, . . .. Thus, p{z) can be identified with the eigenvalue density. 

Although the simple case t = 1 has appeared to a large extent in the liter- 
ature, we will now discuss it to give the reader a better feeling of the results 
and to show the relation to other methods to investigate the distribution of 
eigenvalues. 

So from now on let t = 1 and define ao = a, bo = b/2. Henceforth, we 
put a > and 6 > without loss of generality (it is easily seen that p{z) 
for the pairs a, b and a, —b is the same). It will be interesting to note the 
distinction between the cases a > b and a < b. 

Corollary of the theorem Let t = 1, clq — a, bf) — b/2. Under the 
conditions specified above, the eigenvalue density of J{n)/(f{n) as n ^ oo is 
the following: 

1) 0< a/b < 1, 

p{z) = ^I^/^a+b sign z) ^k-^^j^llZau^ ^ ^ [a~b,a+b], and p{z) = otherwise; 

2) a/b > 1, 



5 



j_ rl g(w)aLU 



g{ijj)du) 



0, 



z & [a — h,a + h] 
[0, a + 6] 



Note that for ^^(n) = ra'^, 7 > we can write case 2 in a simpler form: 

-1+1/7 




z e [0, a - 6] 
2; G [a — 6, a + 6] 
z ^ [0, a + 6] 



where 



7r7 7^/(a+6) y'62^2_ (^_a^)2' 

The form of p{z) is illustrated for ^{n) = rC , in Figures 1,2. When 
7 — i> the density p{z) for all a, 6 approaches the distribution po(2^) = 



,z — a) 



z & [a — b,a + b]. 



One can note that for ^p{n) = n the matrix J/r, where r = y |a^ — 6^1, 
has the same asymptotics of the matrix elements as the matrix associated 
with the Meixner (for a > 6) or Meixner-Pollaczek (for a < b) orthogonal 
polynomials.^ Asymptotics of the eigenvalues of J{n)/r is the asymptotics 
of the zeros (to the main order) of these polynomials. 

For (p{n) = n the function p{z) is elementary. Integration gives in cases 
1 and 2, respectively 



1) P(^) 



2) p{z) 



— arccosh 4: + f 
0, 



-arccos (-^ + f y , 




z E [a — b,a + b] 
z ^ [a — b, a + b] ' 

z e [0, a - 6] 

2; G [a — b,a + b] . 

z ^ [0,a + b] 



< arccosfa;) < vr. 



Here we disregard a sign in front of the matrix and polynomials. 



6 



a-b 



a+b 



Figure 1: Example of the density of eigenvalues p{z) for a < b, (p{n) = rC . 
Solid curve corresponds to some 7 > 1, dashed curve, to < 7 < 1. 

Replacing r by 1 in these formulas gives the contracted asymptotic densities 
of zeros of the Meixner-Pollaczek (case 1) and Meixner (case 2) polynomials. 
The quantities a, h are expressed in this case in terms of the parameters of 
the polynomials. For detailed discussion of zero densities for various systems 



of orthogonal polynomials see and references therein. 

The Meixner and Meixner-Pollaczek polynomials satisfy difference equa- 
tions of the type |jl3| : 

B{x)pn{x + 6) - C{x, ■n)pn{x) + D{x)pn{x - 5) = 0. (11) 

By analyzing such equations, it is possible to obtain information on the 



asymptotic distribution of zeros of Pn{x) |jTj]. Thus, corollary of the theorem 
for ip{n) = n can be given an alternative proof, and also more precise informa- 
tion (than the density) about asymptotics of zeros can be found: distribution 
of zeros of the Meixner polynomials pn{nz) in the region 2; e [0, a — 6] where 
p{z) = 1. 
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a-b 



a+b 



Figure 2: Example of the density of eigenvalues p{z) for a > b, <f{n) = rC . 
Solid curve corresponds to some 7 > 1, dashed curve, to < 7 < 1 (in fact, 
it is obvious from the picture that the second derivative of the dashed curve 
is positive for 2; e [0, a — 6], so it corresponds to some 7 < 1/2). 
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